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ABSTRACT 

In this paper we analyze the stellar populations present in M31 using nine sets of adjacent 
HST-NICMOS Camera 1 and 2 fields with galactocentric distances ranging from 2' to 20'. These 
infrared observations provide some of the highest spatial resolution measurements of M31 to 
date; our data place tight constraints on the maximum luminosities of stars in the bulge of M31. 
The tip of the red giant branch is clearly visible at M\, Q i ~ —3.8, and the tip of the asymptotic 
giant branch (AGB) extends to Aff, ; ~ —5. This AGB peak luminosity is significantly fainter 
than previously claimed; through direct comparisons and simulations we show that previous 
measurements were affected by image blending. We do observe field-to-field variations in the 
luminosity functions, but simulations show that these differences can be produced by blending in 
the higher surface brightness fields. We conclude that the red giant branch of the bulge of M31 
is not measurably different from that of the Milky Way's bulge. We also find an unusually high 
number of bright blueish stars (7.3 arcmin^ 2 ) which appear to be Galactic foreground stars. 

Subject headings: galaxies: individual(M31/NGC224) 
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1. Introduction 

The first deep infrared (IR) observations of 
stars in the bulge of our Galaxy were carried out 
by Frogel & Whitford (1987). Measuring the lu- 
minosity function (LF) using the M giant grism 
surveys of Blanco, McCarthy & Blanco (1984) and 
Blanco (1986), they found many luminous giants, 
but noticed that the LF has a sharp break at 
M bot ~ -4.5 (M K ~ -7.5), with the brightest 
stars extending to M^oi — —5. 

The tip of red giant branch (RGB), defined 
by the core mass required for helium flash, oc- 
curs at a luminosity of Mboi — —3.8. Any 
stars brighter than this limit are therefore on 
the asymptotic giant branch (AGB). The stars 
observed in the Galactic bulge extend ~ 1.2 mag- 
nitudes brighter than the tip of the RGB. Since 
metal-poor ([Fe/H]< —1) Galactic globular clus- 
ters do not exhibit such luminous AGB stars, this 
might have suggested a younger age for the bulge 
population, since the luminosity of the brightest 
AGB stars increases with decreasing age (e.g. Iben 
& Rcnzini 1983). However metal-rich globular 
clusters do have stars which can reach luminosi- 
ties of Mboi — —5.0, while still having ages com- 
parable to the metal poor clusters (Frogel & Elias 
1988; Guarnieri, Renzini & Ortolani 1997). More- 
over, it has been demonstrated that the stellar 
population of the Galactic bulge is dominated by 
metal-rich stars ([Fe/H]> -1) (McWilliam & Rich 
1994) that are as old as Galactic globular clusters 
(Ortolani ct al. 1995; Fcltzing & Gilmore 2000; 
Kuijken & Rich 2002; Zoccali et al. 2002), and 
that the number of stars brighter than the RGB 
tip is consistent with the frequency observed in 
old, metal rich globular clusters. In summary, the 
stellar population in the bulge of the Milky Way 
is as old as the oldest Galactic globular clusters, 
and old, metal rich stellar populations are able to 
produce AGB stars only as bright as M^oi — —5. 

Work was already underway studying stars in 
other nearby galaxies to determine whether the 
properties of the stars in the Galactic bulge are 
typical of all bulges. The nearest and brightest 
large spiral, M31 (the Andromeda Galaxy) was 
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the obvious first choice for a comparison. Some of 
the first measurements of stars in the inner bulge 
of M31 (~ 1 kpc from the nucleus) were made 
by Mould (1986). He found that M31's brightest 
bulge stars were ~ 1 magnitude more luminous 
than the brightest stars in M31's halo. Rich et al. 
(1989) then took spectra of some of these bright 
stars, and found most to have properties charac- 
teristic of late- type M giants. 

It was unclear what was causing the appar- 
ent difference between the stellar populations of 
the bulges of the Milky Way and M31. The de- 
pendence of the AGB peak luminosity on mass 
(and therefore age and mass loss) and metallicity 
pointed to several possible explanations for this 
observed difference. However, all of these expla- 
nations implied a difference in the formation or 
evolutionary processes of these two otherwise very 
similar galaxies. 

To see if the luminous stars in M31 are in- 
deed similar to those found in the bulge of the 
Milky Way, Rich & Mould (1991) measured a sam- 
ple of ~ 600 stars in the inner bulge of M31 
4' from the nucleus with the Hale 5m reflector 
(see §8.2). Their resulting LF had a drop at 
M}, i ~ —4.5, similar to that seen in BW, but ex- 
tended to Mboi ~ —5.5. To explain these excess 
luminous stars they proposed several theories. (1) 
These stars could be younger stars from the disk 
superposed on the bulge. (2) They could be super- 
metal-rich in chemical composition, since the lu- 
minosity of the brightest AGB stars increases with 
metallicity. (3) The bulge of M31 could have a 
young stellar component, since the luminosity of 
the brightest AGB stars increases with decreasing 
age (Iben & Renzini 1983). (4) They could be the 
result of merged lower-mass main-sequence stars 
(blue straggler progenitors), which produces more 
massive and luminous stars. 

Soon thereafter, Davies, Frogel & Terndrup 
(1991) imaged the M31 bulge in the near-IR us- 
ing the 3.8m UKIRT facility, albeit 7.2' from the 
nucleus, a factor of two more distant than Rich. 
Their LF has an upper limit ~ 0.5 mag brighter 
than that seen in the Galactic bulge. They argued 
that contamination by stars from a young disk, 
which lies behind the bulge, is most likely respon- 
sible for the luminous stars observed by RM91. 

DePoy et al. (1993) carried out a if-band sur- 
vey of ~ 17000 stars in 604 arcmin 2 of Baade's 
Window to check the possibility that the M gi- 
ant surveys in the bulge of the MW (Blanco, Mc- 
Carthy & Blanco 1984; Blanco 1986) may have 
missed very luminous stars similar to the ones seen 
in the bulge of M31. The DePoy et al. (1993) 
observations turned up no such population of lu- 
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minous stars, and their derived LF is consistent 
with that obtained by Frogel & Whitford (1987). 
In order to compare their Galactic Bulge obser- 
vations with those of the bulge of M31, DePoy et 
al. (1993) rebinned and smoothed their image to 
simulate the M31 observations. The resulting de- 
graded image showed that few if any of the "stars" 
on this simulated M31 field corresponded to indi- 
vidual stars on the original image; most were just 
random groupings of stars. A quantitative anal- 
ysis showed that the extreme crowding caused an 
artificial brightening in the LF of more than one 
magnitude. They thus concluded that the lumi- 
nous stars seen in M31's bulge were most likely 
not real, but an artifact of image crowding. 

At nearly the same time Rich, Mould & Gra- 
ham (1993) acquired new observations of five fields 
in the bulge of M31 with the Palomar IR Imager 
on the Hale 5m telescope (see §8.1). By measur- 
ing the LFs in fields with different expected disk 
contributions from 2 to 11 arcminutes from the 
nucleus, they rejected the hypothesis of Davies, 
Frogel & Terndrup (1991) that the bright stars 
are disk contaminants. Using their own model 
fields, which showed that they could accurately 
measure the GB tip despite the crowding, they 
also argued against the idea of DePoy et al. (1993) 
that the bright stars are stellar blends. While 
they did concede that some of their measurements 
may have been affected by crowding of up to 1 
magnitude, they maintained that they were not 
generally measuring clusters of blended images. 
Further calculations by Rich, Mould & Graham 
(1993) also showed the numbers of blue straggler 
progeny stars to be insufficient to explain the num- 
ber of luminous stars. 

In pursuit of a resolution to this controversy, 
Rich & Mighell (1995) obtained HST Wide-Field 
Planetary Camera (WFPC1) observations of the 
inner bulge of M31. Regretfully, these observa- 
tions were taken with HST's original aberrated 
optics, reducing the effective resolution to barely 
better than was available from the ground. These 
observations also yielded many luminous stars, al- 
though not quite as bright as previously measured. 
The data also suggested that the brightest stars 
may be concentrated toward the center of M31. 

Approaching this problem from the theoretical 
side, Renzini (1993, 1998) performed calculations 
to estimate the number of stars in all evolutionary 
stages in each pixel. He showed that the num- 
ber of blends increases quadratically with both 
the surface brightness of the target and with the 
angular resolution of the observations. Applying 
these calculations to existing photometric data for 
the inner bulge of M31, he concluded that all pre- 



vious ground-based observations were dominated 
by blends, and even questioned the HST obser- 
vations of Rich & Mighell (1995), pointing to the 
measured blue (R — I) colors as indicative of their 
blended origin. 

With the Wide-Field Planetary Camera 2 
(WFPC2), Jablonka et al. (1999) observed three 
fields in the bulge of M31 at optical wavelengths. 
With HST's improved resolution, they did not find 
stars more luminous than those in the Galactic 
Bulge, and concluded that previous detections of 
very bright stars were likely the result of blended 
stars due to the crowding in WFPC1 and ground- 
based images. 

However, since even very luminous evolved 
stars can go undetected at optical wavelengths 
due to molecular blanketing, Davidge (2001) re- 
cently obtained new infrared images of the bulge 
of M31 with the 3.6m CFHT. With the help of 
adaptive optics (AO), his JHK observations con- 
firm the optical non-detection of very luminous 
stars made by Jablonka et al. (1999). Although 
there is agreement between the brightest stars 
measured by (Davidge 2001) in the bulge of M31 
and the brightest stars measured in the Galactic 
Bulge (Frogel & Whitford 1987), the luminosity 
functions still show considerable differences. The 
M31 bulge LF measured by Davidge (2001, Fig- 
ure 7) does not show a break at Mk ~ —7.5 as 
is observed in BW, but instead shows a change in 
slope at Mk ~ —8 and a break at Mk ~ —8.6, 
indicative of different star formation histories in 
the MW and M31 if correct. 

Thus it appears that this decade-old con- 
troversy has not yet been completely resolved. 
While there now appears to be agreement that 
the brightest of the previous measurements were 
blends, it is still not certain whether the lumi- 
nosity functions of the bulges of the Milky Way 
and M31 are consistent with one another. In this 
paper we will show that indeed they are, and that 
even the most recent observations, including our 
own, are still affected by blending in the inner 
regions of M31. 

The layout of the paper is as follows. We start 
by describing our observations in Section 2, and 
our reduction and photometric techniques in Sec- 
tion 3. We present the color magnitude diagrams 
and luminosity functions in Sections 4 and 5 re- 
spectively. Section 6 gives a brief theoretical anal- 
ysis of blending in M31, followed by detailed sim- 
ulations of all our fields in Section 7. In Section 8 
we compare our measurements with previous ob- 
servations; Rich, Mould & Graham (1993) in §8.1, 
Rich & Mould (1991) in §8.2, and Davidge (2001) 
in §8.3. Section 9 discusses the bright stars which 



3 



appear to be Galactic foreground stars, and we 
conclude with a brief summary of our results in 
Section 10. 

2. Observations 

In this paper we analyze images of M31 taken 
from two different NICMOS proposals. Proposal 
7876 imaged the five central fields of Rich, Mould 
& Graham (1993). These fields were carefully cho- 
sen to sample varying bulge/disk ratios, and are 
indicated by F1-F5 on Figure 1. Proposal 7826 
imaged five globular clusters in M31. The four 
fields used in this study are indicated on Figure 1 
by their cluster numbers: F170, F174, F177 and 
F280. We omit the Gl field from this analysis, 
since at 34 kpc from the center of M31, the frame 
is dominated by cluster stars. 

Since the NIC1 and NIC2 cameras are at dif- 
ferent positions in the HST focal plane, we can si- 
multaneously image 2 fields at each pointing (see 
Figure 2). The NIC2 field is 19.2" across, and sep- 
arated by 17.5" from the 11" NIC1 field (32.6" be- 
tween field centers). Their different sizes are due 
to their different spatial resolutions; the NIC1 im- 
ages have a plate scale of 0.043" pixel -1 , compared 
to the 0.0757" pixel -1 of NIC2. Having different 
resolutions at each pointing will prove useful for 
understanding the severe crowding very near the 
center of M31. 




RA (deg) 



Fig. 2. — Map of our observations. The filled 
squares represent our NIC1 (small) and NIC2 
(large) fields. The large open squares show the ob- 
servations of Rich, Mould & Graham (1993). The 
contours are r-band surface brightness contours 
taken from Kent (1989) and interpolated to 0.5 
magnitude intervals. This figure illustrates the rel- 
ative positions, sizes and separations of the NIC1 
and NIC2 fields, as well as why it is not completely 
fair to compare parallel NIC1 and NIC2 fields un- 
der the assumption that they have equal surface 
brightness. 

Our observations are summarized in Table 1. 
The top half of the table lists our NIC1 observa- 



tions, and the second half is for NIC2. The first 
column lists the field ID; the second and third give 
the field coordinates. The distance from the cen- 
ter of M31 in arcminutes is listed in column 4. 
Using this distance, the position angle from the 
major axis of M31, the r-band surface brightness 
from Kent (1989), and an assumed (r — K) color 
of 2.9, we estimate the i-T-band surface brightness 
of each field, which is listed in column 5. Taking 
Kent's bulge - disk decomposition, we also give 
the bulge/disk ratio in column 6. 

The NICMOS focus was set at the compromise 
position 1-2. This is the best focus for simultane- 
ous observations with cameras 1 and 2. All of our 
observations used the MULTIACCUM mode (MacK- 
enty et al. 1997) because of its optimization of the 
detector's dynamic range and cosmic ray rejection. 

Our primary (pointed) observations were taken 
with NIC2. Each field was observed through 
three filters: F110W (0.8-1.4 fj,m), F160W (1.4- 
1.8 ,uni), and F222M (2.15-2.30 /an). These filters 
are close to the standard ground-based J, H, & K 
filters. However, to maximize the depth of our par- 
allel NIC1 observations, we only used the F110W 
(J) filter. Total integration times and FWHMs 
for each camera and filter combination are given 
in Table 2. The dates of the observations are listed 
in Table 3. 

All the observations implemented a spiral 
dither pattern with 4 positions to compensate 
for imperfections in the infrared array. For fields 
F1-F5 the dither steps were 0.4" for the J and 
K band images, and 1.0" for the H band images. 
Fields F170-F280 are the same, except that we 
used 5.0" dithers in H . 

As previously mentioned, the observations of 
fields F170-F280 are from another proposal tar- 
geting M31's metal-rich globular clusters (Stephens 
et al. 2001b). In these observations we exclude 
stars inside radii of 1.4", 0.5", 0.6", and 5.0" 
around the clusters G170, G174, G177 and G280 
respectively, to avoid cluster stars. 

Images of fields F1-F5 are shown in Figures 3 
and 4 for NIC2 and NIC1 respectively. The NIC1 
images of fields F170-F280 are shown in Figure 5, 
and their NIC2 counterparts are given in Figure 
2 of Stephens et al. (2001b). These images are 
the combination of 4 and 12 dithers for NIC2 and 
NIC1 respectively. The dimensions of each set of 
combined images are different due to the varying 
plate scale and dither size, and are given in the 
figure captions. 

When converting to absolute or bolometric 
magnitudes we assume a distance modulus to 
M31 of (m - M) = 24.4, which corresponds 
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Table 1 
M31 Observations 



ID 




<5 a 


Radius 




Rill ova /~ni«lr 


y- 1 -) 


(21 


(3) 


(41 


(5) 


(6) 


NIC1 












Fl 


h 42 m 37!83 


41°14'18.3" 


2.20' 


15.09 


6.4 


F2 


h 43 m 28!46 


41°6'48.9" 


12.49' 


18.72 


0.1 


F3 


h 42 m 32!63 


41°12'45.8" 


4.03' 


15.91 


3.5 


F4 


h 43 m 2U3 


41°13'2.0" 


4.57' 


16.88 


1.3 


F5 


h 42 m 25!51 


41°11'39.0" 


5.72' 


16.37 


2.1 


F170 


h 42 m 32!14 


41°9'51.5" 


6.69' 


16.68 


1.7 


F174 


h 42 m 33U5 


41°17'4.5" 


2.29' 


15.52 


4.9 


F177 


h 42 m 35?30 


41°13'27.4" 


3.17' 


15.58 


4.7 


F280 


h 44 m 30?63 


41°21'0.7" 


20.55' 


18.35 


0.2 


NIC2 












Fl 


h 42 m 36!30 


41°14'51.7" 


1.97' 


14.98 


7.7 


F2 


h 43 m 27?07 


41°7'23.0" 


11.89' 


18.54 


0.1 


F3 


h 42 m 30?87 


41°13'17.7" 


3.80' 


15.82 


3.9 


F4 


h 43 m 0?44 


41°13'34.3" 


3.98' 


16.57 


1.9 


F5 


h 42 m 22!63 


41°11'57.8" 


5.84' 


16.40 


2.0 


F170 


h 42 m 32?40 


41°10'29.0" 


6.08' 


16.52 


1.9 


F174 


h 42 m 33?30 


41° 17' 42.0" 


2.59' 


15.78 


4.1 


F177 


h 42 m 34?40 


41°14'3.6" 


2.79' 


15.40 


5.5 


F280 


h 44 m 29?50 


41°21'36.0" 


20.49' 


18.33 


0.2 



a J2000 

b mag/arcsec~ 2 , from Kent (1989) assuming (r — K) = 2.9. 



Table 2 

NICMOS Exposure TiMES a & FWHMs 





NIC1 




NIC2 




Fields 


F110W 


F110W 


F160W 


F222M 


F1-F5 


4992 


1280 


2048 


1664 


F170-F280 


7552 


1920 


3328 


2304 


FWHM 


0.099" 


0.125" 


0.148" 


0.185" 



a Exposurc times in seconds. 



Table 3 
Observation Dates 



Field 


Date 


Fl 


1998/09/20 


F2 


1998/09/18 


F3 


1998/09/24 


F4 


1998/09/23 


F5 


1998/10/13 


F170 


1998/08/10 


F174 


1998/08/13 


F177 


1998/09/08 


F280 


1998/09/13 
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Fig. 1. — A 1° field from the Digitized Sky Survey showing the location of the nine fields used in our analysis. 
North is up and East is to the left. 



to 3.7 parsecs arcsec and an extinction of 
E(B — V) = 0.22. 

3. Data Reduction & Photometry 

Our data were reduced with the STScI pipeline 
supplemented by the Image Reduction and Anal- 
ysis Facility (IRAF 6 ) nicproto package (May 
1999) to eliminate any residual bias (the "pedestal" 
effect). Object detection was performed on a com- 
bined image made up of all the dithers of all the 
bands (12 images in total). Point spread func- 
tions (PSFs) were determined from each of the 
four dithers, then averaged together to create 
a single PSF for each band of each target. In- 
strumental magnitudes were measured using the 
ALLFRAME PSF fitting software package (Stetson 
1994), which simultaneously fits PSFs to all stars 
on all dithers. DAOGROW (Stetson 1990) was then 
used to determine the best magnitude in a 0.5" 
radius aperture. 

6 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by AURA, Inc., under 
cooperative agreement with the NSF. 



We finally transformed our photometry to the 
CIT/CTIO system. The NIC2 measurements used 
the transformation equations of Stephens ct al. 
(2000), listed in equations 1-3. The NIC1 trans- 
formation proved to be more complicated, and is 
based on a comparison with the much lower spa- 
tial resolution groundbased observations of Rich, 
Mould & Graham (1993). A detailed description 
of the technique is given in Appendix A. 

The selection criteria are also slightly different 
for stars measured in the NIC1 and NIC2 fields. 
For the NIC2 frames we require measurements in 
all three bands, with PSF-fitting errors smaller 
than 0.25 magnitudes in each band. For NIC1, 
with only J-band observations, we only require 
that the PSF-fitting error be less than 0.25 mag- 
nitudes. 



mj = TOno - (0.198 ± 0.036)(to 110 - m 2 2 2 ) 

+ (21.754 ±0.030) (1) 

m H = mien - (0.177 ± 0.037)(toiio - m 22 2) 

+ (21.450 ±0.028) (2) 

m K = 777 222 + (0.074 ± 0.037)(toiio - m 222 ) 




. 3. — NIC2 J-band images of Fields F1-F5, each image is ~ 20" across. 
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4. — NIC1 J-band images of Fields F1-F5, each image is ~ 12" across. 
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5. — NIC1 J-band images of Fields F170-F280, each image is ~ 16" across. 
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+ (20.115 ±0.031) (3) 
4. Color-Magnitude Diagrams 

The Mko~(J — K)o color-magnitude diagrams 
of all 9 NIC2 fields are shown in Figure 6. The 
overplotted solid lines are contours of constant 
bolometric magnitudes of —4 and —5, based on 
the bolometric corrections for M giants in Baade's 
Window calculated by Frogel & Whitford (1987). 
These plots assume a distance modulus to M31 
of (to - M) = 24.4, E(J - K) = 0.12, and 
A K = 0.07. 

The RGB and AGB are both visible in these 
CMDs. The tip of the RGB is more clearly defined 
in the less crowded fields, and a differential bolo- 
metric LF shows that it occurs at M^oi ~ —3.75. 
In the more crowded fields the RGB tip gets 
blurred because of blending, which pushes stars 
up off of the RGB. The tip of the bulge AGB ex- 
tends to Mboi ~ — 5. This AGB tip is significantly 
fainter than previously claimed, and we address 
this in a comparison with previous observations in 
Section 8. 

4.1. LPVs 

Long-period variables (LPVs) are large ampli- 
tude, luminous red variable stars with periods 
ranging from 50 to several hundred days. These 
stars are on the AGB and represent the brief fi- 
nal stages of low- to intermediate-mass stellar evo- 
lution on the giant branch. Based on measure- 
ments of variables in the Galactic Bulge (Frogel 
& Whitford 1987), we have marked the region 
of each CMD where we expect to find primar- 
ily LPVs (Fig. 6). This region is indicated by 
a dashed box in the upper right of each CMD, 
with (J — K)o > 1.6 and Mko < —6. Since im- 
age blending generally shifts objects to bluer col- 
ors (Stephens et al. 2001a), this region should be 
relatively insensitive to crowding causing spurious 
LPV candidates, although in extreme cases, LPVs 
may actually be shifted blucward out of the box, 
thus only giving us a lower limit to the number 
of LPVs. A casual comparison of this LPV re- 
gion between fields shows that some of the CMDs, 
particularly from the inner fields, have many more 
potential LPVs. 

Since the relative numbers and luminosities of 
LPVs are sensitive to the age and/or metallicity 
of the parent population, LPVs can be used to 
look for field-to-field variations in the stellar pop- 
ulation. Frogel & Whitelock (1998) argue that 
the relative number of LPVs to non-variable gi- 
ants is independent of [Fe/H] for [Fe/H]< 0, and 
for higher metallicities the LPV lifetime is signif- 
icantly reduced due to increased mass loss rates. 



Thus for stellar systems with a super-solar metal- 
licity component, the ratio of non- variables to 
LPVs will appear high compared to lower metal- 
licity systems. This idea seems to be supported 
by their determination of a higher non- variable gi- 
ants to LPV ratio in the Galactic bulge compared 
to globular clusters. 

To determine whether there is a change in the 
relative numbers of LPV candidates among our 
fields, we compare their numbers to the number of 
nonvariablc giants, classified by having (J — K)q < 
1.6 and Mko < — 5. The ratio of giants to LPVs 
for each of our NIC2 fields is shown in Figure 7 
as a function of bulge/disk ratio. This plot shows 
no trend in the giant/LPV ratio, instead all of our 
observations are scattered around the average ra- 
tio of N(Giants)/N(LPVs)~ 24. Thus the LPVs 
appear to be uniformly distributed with the non- 
variable stars. 

100 I 

80 - 



60 



20 



10 



Bulge / Disk 

Fig. 7. — The ratio of the number of nonvariablc 
giants to long-period variable (LPV) candidates as 
a function of bulge-to-disk ratio. Classifications 
as an LPV or non-variable giant are based on the 
stellar colors and luminosities indicated on each 
CMD in Figure 6. All fields are consistent with 
the average ratio of 24. The errorbar with the 
lowest bulge/disk ratio belongs to field F2, which 
has 80 giants, but no LPVs. 
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Fig. 6. — The dereddened color magnitude diagrams for each of the 9 NIC2 fields. We have drawn on lines 
of constant bolometric magnitude at M^oi — —4 and —5, using the bolometric corrections calculated for 
Baade's Window M giants by Frogel & Whitford (1987). The box in the upper right of each panel, with 
(J — K) > 1.6 and Mko < —6, indicates the region we expect to find primarily LPVs. The box in the upper 
left of each panel, with (J — K) < 1.6, and Mko < — 5, is the region we use to count nonvariable giants. 
We have assumed E(B -V) = 0.22, giving E(J - K) = 0.12 and A K = 0.07. 
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5. Luminosity Functions 

The J-band LFs measured in all 18 fields are 
shown in Figure 8 and listed in Table 4 in units 
of number per magnitude per square arcsecond. 
The figure shows both the NIC1 and NIC2 LFs 
overplotted for each field. The first thing to note 
about this compilation of M31 LFs is that the 
NIC1 and NIC2 measurements are not exactly the 
same. While most show good agreement, several 
of the NIC2 LFs extend to brighter magnitudes, 
and all of the NIC1 LFs extend fainter. 

The faint end differences are a combined re- 
sult of the longer exposure times and better spa- 
tial resolution of NIC1. The NIC1 exposure times 
are nearly four times longer than those of the J- 
band NIC2 observations. This is because while 
NIC2 was cycling through all three J,H & K fil- 
ters, NIC1 observed only through the J-band fil- 
ter. The faint end photometry is also affected by 
the level of crowding in the field. NIC2 is under- 
sampled at J, making it more difficult to distin- 
guish close objects. Thus in very crowded fields, 
NIC1 has an advantage over NIC2, accounting for 
the larger faint-end difference seen in the more 
crowded fields. 

There are several reasons for the differences 
seen at the bright end. The first is the difference 
in the field of view of each camera. NIC2 cov- 
ers an area on the sky three times that of NIC1, 
and thus has a much better chance of finding the 
rarer brighter stars. Second is resolution and im- 
age sampling. NIC1 has a well sampled J-band 
PSF with a FWHM of 0.099", while the NIC2 J- 
band PSF is undcrsampled with a larger FWHM of 
0.125". Also, the NIC2 photometric calibration re- 
quires detections in all three bands. Thus the lim- 
iting resolution for NIC2 is actually the K-b&nd, 
which has a FWHM of 0.185", nearly twice the size 
of the NIC1 J-band PSF. The combined result is 
that the NIC2 observations are more sensitive to 
blending, which can artificially brighten stars, ex- 
tending the bright end of the LF. Finally, there 
is the coincidence that most of the NIC1 obser- 
vations occur in regions of slightly fainter surface 
brightnesses than their NIC2 counterparts (see Ta- 
ble 1 and Figure 2), somewhat exacerbating both 
aforementioned effects. 

The JT-band NIC2 LFs are shown in Figure 
9 and listed in Table 5, both normalized to give 
number/arcsec 2 /magnitude. Although NIC1 is 
more resilient against blending, both NIC1 and 
NIC2 observations are susceptible to its ill effects 
(see §7). However, under the hypothesis that all 
the measured bulge LFs arise from a single true 
LF, and that the differences are purely a result 



of blending, we estimate that the tip of the AGB 
occurs at Mj ~ —6.25 and Mk ~ —8. 

5.1. Comparison with Baade's Window 

We show our NIC2 if -band luminosity func- 
tions superposed on the LF of the Galactic Bulge 
in Figure 10. The Bulge LF is a composite of mea- 
surements made in Baade's Window (BW), with 
the bright end (M K < —6.6) from Frogel & Whit- 
ford (1987) and the faint end from DePoy et al. 
(1993). All the LFs have been normalized in the 
range -7 < M K < -5.5. 

The bottom panel of Figure 10 is a compari- 
son between BW and the seven bulge fields in our 
sample. These fields all have bulge-to-disk ratios 
greater than one, as listed in Table 1. These fields 
are all of very high surface brightness, and there- 
fore we expect most to exhibit some amount of ar- 
tificial brightening due to blending. However, as 
the plot shows, there is still very good agreement 
between these M31 bulge fields and the bulge of 
the MW. 

Thus based on our infrared luminosity func- 
tions, the stellar population of bulge of M31 is very 
similar to that of the Milky Way. The match up 
between the LFs is quite good, and when one takes 
into account our prediction of a small amount of 
artificial brightening due to blending in our most 
crowded bulge fields, the correspondence will be 
even better. 

5.2. Comparison of Disk and Bulge LFs 

The top panel of Figure 10 is a comparison be- 
tween BW and the two disk fields in our sample. 
The F2 & F280 fields both have bulge-to-disk ra- 
tios less than one (see Table 1). These two fields 
also have the lowest surface brightnesses, which 
means that the effects of blending are the least, 
and hence their photometric measurements are the 
most trustworthy. However in this comparison, 
we see that both of these disk LFs extend slightly 
brighter than the break measured in BW, and do 
so more prominently than any of the bulge fields. 

In order to determine whether or not the mea- 
sured disk and bulge luminosity functions arc in 
fact distinguishable from one another, we com- 
pare the distributions of stellar luminosities using 
the KS-tcst. We combine the F2 and F280 mea- 
surements to represent the disk population, and 
use all other fields for the bulge population. In- 
completeness in the more crowded fields limits the 
comparison to only bright stars with < —6, 
and in this range the KS-test shows a conspicu- 
ous overabundance of luminous disk stars in the 
range —8 < Mk < —7.5. However the signifi- 
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Table 4 
J- band Luminosity Functions 



NICl NIC2 



J 


Fl 


F2 


F3 


F4 


F5 


F170 


F174 


F177 


F280 


Fl 


F2 


F3 


F4 


F5 


F170 


F174 


F177 


F280 


15.875 


0. 


.00 


0. 


00 


0. 


00 





.00 





.00 





.00 


0.00 


0.00 





.00 





.00 


0. 


.00 





.00 





.01 





.00 





.00 





.00 


0. 


00 


0.00 


16.125 





00 





00 





00 





.00 





.00 





.00 


0.00 


0.00 





.00 





.00 





00 





.00 





.00 





,00 





.00 





00 





00 


0.00 


16.375 





.00 


0. 


00 


0. 


00 





.00 





.00 





.02 


0.00 


0.00 





.00 





.00 





.00 





.00 





.00 





.00 





,00 





.00 





00 


0.00 


16.625 


0. 


00 





00 


0. 


00 





.00 





.00 


0. 


.00 


0.00 


0.00 





.00 





.00 





00 





.00 





.01 





.01 





.00 


0. 


00 


0. 


00 


0.00 


16.875 





00 


0. 


.00 


0. 


00 





.00 





.00 


0. 


.00 


0.00 


0.00 





.00 





.00 





00 





.01 





.00 





.00 





.00 


0, 


00 


0. 


00 


0.00 


17.125 





.00 


0. 


00 





00 





.00 





.00 





02 


0.00 


0.00 





.00 





.01 





.01 





.00 





.00 





.00 





.00 


0. 


.00 


0. 


00 


0.00 


17.375 





00 





00 


0. 


00 





.00 





.00 


0. 


.00 


0.00 


0.02 





.00 





.02 





01 





.01 





.00 





.00 





.00 


0. 


01 


0. 


00 


0.00 


17.625 





.03 


0. 


.00 


0. 


00 





.00 





.00 


0. 


.00 


0.00 


0.02 





.00 





.02 





.00 





.02 





.00 





.00 





.00 


0. 


.03 


0. 


03 


0.01 


17.875 





11 


0. 


00 


0. 


00 





.00 





.00 


0. 


.02 


0.02 


0.06 





.02 





16 





00 





.02 





.02 





.05 





.02 


0. 


09 


0. 


12 


0.00 


18.125 


0. 


11 





03 


0. 


06 





.00 





.03 


0. 


.03 


0.17 


0.16 





.05 





.48 





01 





15 





.02 





.06 





.06 


0. 


08 


0, 


17 


0.01 


18.375 





.19 


0. 


.00 





06 





.11 





.03 


0. 


.08 


0.20 


0.12 





.02 


1 


.01 





.01 





17 





.07 





.08 





,11 


0. 


.18 





52 


0.06 


18.625 


0. 


92 


0. 


.03 


0. 


.19 





.00 





.17 


0. 


.09 


0.39 


0.39 





.02 


1 


.70 


0. 


.03 





.38 





.08 





.25 





,11 


0. 


.56 


0. 


76 


0.03 


18.875 


2. 


.36 





00 





22 





.36 





.14 





.22 


0.89 


0.70 





.02 


3 


.95 





.02 





.92 





.39 





.50 





,37 


1, 


.34 


2. 


00 


0.01 


19.125 


4. 


14 


0. 


11 


1. 


11 





36 


1 


.08 


0. 


.59 


2.59 


2.31 





.05 


4 


.90 





08 


1 


.79 





.81 


1 


.04 





,81 


2, 


28 


3, 


01 


0.12 


19.375 


5. 


.81 


0. 


03 


2. 


33 





.72 





.94 


0. 


91 


3.14 


3.30 





.14 


5 


33 


0. 


.11 


2 


.39 


1 


.04 


1 


.32 


1 


,01 


2. 


.71 


3 


45 


0.09 


19.625 


6 


.00 


0. 


22 


2. 


.72 


1 


.47 


1 


.72 


1 


.05 


4.14 


3.86 





.09 


5 


.12 





.08 


2 


.83 


1 


.41 


1 


.50 


1 


19 


2. 


.86 


4. 


13 


0.06 


19.875 


7. 


08 





14 


2. 


.89 


1 


.22 


1 


.58 


1. 


27 


3.67 


4.16 





.14 


5 


.31 


0. 


19 


2 


.67 


1 


.27 


1 


.42 


1 


,57 


2. 


76 


3, 


80 


0.13 


20.125 


7. 


.69 


0. 


25 


3. 


39 


1 


.22 


1 


.92 


0. 


.98 


4.33 


4.89 





.19 


5 


.24 


0. 


.24 


3 


.00 


1 


.60 


1 


.69 


1 


.25 


3, 


.23 


4, 


19 


0.10 


20.375 


9 


.33 


0. 


03 


3. 


61 


1 


.42 


1 


.50 


1 


60 


5.52 


5.47 





.23 


4 


.44 





.18 


3 


.59 


1 


.68 


1 


.87 


1 


.66 


3, 


48 


4, 


05 


0.16 


20.625 


10.28 





14 


4 


86 


1 


.75 


2 


.33 


1. 


61 


5.91 


6.00 





.17 


3 


.61 


0. 


17 


4.03 


2 


.14 


2 


.29 


1 


.69 


3 


48 


3, 


90 


0.19 


20.875 


11.36 


0. 


19 


4. 


.78 


2. 


.06 


3 


,31 


2 


33 


7.69 


6.67 





.34 


2 


.25 


0. 


33 


4 


.35 


2 


.41 


3 


.10 


2 


.30 


3. 


.78 


3 


31 


0.27 


21.125 


12.50 


0. 


36 


6 


25 


2. 


.17 


3 


.97 


2. 


.56 


8.77 


8.34 





33 


1 


.19 





34 


3 


.97 


2 


.68 


3 


.12 


2 


.65 


2. 


90 


1, 


76 


0.41 


21.375 


12.53 


0. 


44 


7. 


61 


3 


.17 


4 


.58 


3 


.33 


8.50 


8.62 





.48 





.51 





37 


2 


.74 


2 


.75 


2 


.79 


2 


,66 


2. 


21 


0, 


88 


0.30 


21.625 


9 


58 





.50 


9 


25 


4 


.06 


5 


.94 


3 


.58 


7.88 


8.36 





.67 





17 


0. 


30 


1 


.60 


2 


.35 


2 


.29 


2 


39 


0, 


70 


0. 


27 


0.45 


21.875 


6. 


36 


0. 


42 


10.56 


5. 


.03 





.83 


4 


.59 


5.41 


5.19 





.56 





.06 


0. 


.35 





.48 


1 


.58 


1 


.44 


1 


,69 


0, 


25 


0, 


05 


0.41 


22.125 


3 


00 


0. 


94 


8.89 


6 


.08 


7 


.47 


3 


.69 


3.12 


3.48 





.64 





.00 





28 





.11 





.68 





.44 





.78 


0, 


06 





01 


0.41 


22.375 


1. 


69 


1. 


00 


6.69 


6 


.14 





.28 


2 


.44 


1.69 


1.97 


1 


12 





.00 





19 





.00 





.19 





.16 





.34 


0. 


01 


0. 


00 


0.20 


22.625 





.47 


1. 


.03 


4.36 


6 


.72 


5 


36 


1 


.28 


0.89 


1.09 





.97 





.00 


0. 


.14 





.00 





.08 





.03 





.09 





.00 


0. 


00 


0.22 


22.875 





.36 


1. 


06 


2.69 


5 


.78 


4 


.28 


0. 


.73 


0.56 


0.36 





.48 





.00 


0. 


07 





.00 





.00 





.00 





.00 


0. 


00 





00 


0.06 


23.125 





19 


1. 


03 


1. 


.86 


4 


.44 


2 


.28 





36 


0.12 


0.20 





.14 





.00 





02 





.00 





.00 





,00 





.00 





00 


0. 


00 


0.01 


23.375 





.03 


0. 


.81 


1. 


11 


3 


.42 


1 


.81 


0. 


.09 


0.11 


0.11 





.09 





.00 





.00 





.00 





.00 





,00 





.00 





00 





00 


0.00 


23.625 





00 


0. 


.58 


0.50 


1 


.92 





.56 


0. 


.06 


0.03 


0.09 





.02 





.00 





00 





.00 





.00 





,00 





.00 


0. 


.00 


0. 


00 


0.00 


23.875 





.00 


0. 


28 


0.08 


1 


.00 





.33 





06 


0.00 


0.00 





.02 





.00 





.00 





.00 





.00 





,00 





.00 


0. 


.00 


0. 


00 


0.00 


24.125 


0. 


.00 


0. 


03 


0.08 





.53 





.06 





.02 


0.00 


0.00 





.00 





.00 





.00 





.00 





.00 





,00 





,00 





.00 


0, 


00 


0.00 




Fig. 8. — The 18 ./-band luminosity functions measured for our 9 sets of NIC1 (beaded) and NIC2 (solid) 
observations. All of the LFs have been normalized to show the number of stars per square arcsec per 
magnitude. The NIC1 LFs go deeper because of slightly better resolution and nearly four times longer 
exposures. We have assumed a distance modulus of (m — M) = 24.4. 



cance is low, with P = 0.34, indicating that the 
two populations are nonetheless consistent with 
being drawn from the same parent population. 

On the other hand, if we limit the comparison 
to only the AGB {Mk < —7), the excess of lumi- 
nous disk stars is enough to drop the KS P— value 
to 0.02. This low probability is marginally signif- 
icant, but is based on a much smaller number of 
stars (12 disk stars and 729 bulge stars). It is also 
noteworthy that the simulations (§7) show no such 
enhancement. 

In summary, both of the two disk fields have a 
slight excess of luminous stars with — 8 < Mk < 
— 7.5, although only statistically significant when 
compared to the bulge fields over a small range 
in luminosity. This small overabundance of AGB 
stars just above the Baade's Window LF break is 
due to the presence of younger disk stars in these 
two fields. 

6. Blending 

In order to analytically estimate the effects of 
blending on our observations, we have used the 
equations of Rcnzini (1998) to predict number of 
stars in each evolutionary stage per resolution ele- 
ment. Several parameters in this calculation have 



a weak dependence on the age, metallicity, and 
IMF of the assumed stellar population. For these 
parameters, we choose to use a ratio of total to K- 
band luminosity of Lt/Lk — 0.36, and a specific 
evolutionary flux B{t) = 2.2 x 10~ n stars yr _1 
L _1 , both suitable for a solar- metallicity, 15 Gyr 
old population. 

To try to estimate the importance of blend- 
ing on fields at different distances from M31, wc 
have calculated the number of RGB stars within 1 
magnitude of the RGB tip, N(RGBT). Since the 
brightest stars in our fields are only ~ 1 magni- 
tude brighter than the expected tip of the RGB, a 
blend with even one RGBT star will distort their 
measurement by > 16%. 

The results of these calculations are displayed 
in Figure 11. The left panel shows the number 
of RGBT stars per resolution element as a func- 
tion of surface brightness for four different imaging 
resolutions. The 0.15" resolution roughly corre- 
sponds to NICMOS, 0.35" is the resolution of the 
Davidge (2001) observations, and 1" corresponds 
to the resolution obtained by Rich, Mould & Gra- 
ham (1993). As an example of using this plot, 
consider an image with 0.35" resolution taken at 
a location where the K-band surface brightness is 
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14.5 magnitudes arcsecond . That image would 
have, on average, 1 RGBT star in each resolution 
element; obviously not conditions favorable for ac- 
curate photometry. 

To make this plot easier to interpret, on the 
right side of Figure 11 we show the number of 
RGBT stars per resolution clement function 
of position in M31, using the same four imaging 
resolutions. To convert from the if-band surface 
brightness in the left panel to radius in the right 
panel, we have used the r-band surface brightness 
measurements of Kent (1989), and we assume that 
(r — K) = 2.9. Following through with the pre- 
vious example, we see that with 0.35" resolution, 
we would find approximately 1 per resolution el- 
ement at a distance of ~ 1.3' from the center of 
M31 along its major axis. 

The question then becomes, what is the limit 
for "good" photometry? While it seems quite ob- 
vious what would be bad, i.e. one or more RGBT 
star per resolution element, what is "good" is more 
difficult to quantify, and requires knowledge of ex- 
actly how good "good" needs to be. Clearly, if 
N(RGBT) is less than one, N(RGBT) is approx- 
imately the probability that a star within one mag- 
nitude of the RGB tip will fall in any given reso- 
lution element. Thus N(RGBT) 2 is the probabil- 
ity that a resolution element will contain a blend 



of two RGBT stars. Therefore the number, and 
severity of the blends which can be accepted de- 
termine the limiting surface brightness. Stephens 
et al. (2001a) have run simulations on their NIC- 
MOS photometry of globular clusters in M31, and 
find that for accurate photometry of stars down 
to K ~ 21, iV(i?GBT)/resolution clement should 
be < 0.05. While this is a good guide, to bet- 
ter understand the observations at different reso- 
lutions and surface brightnesses, it is best to per- 
form simulations to attempt to quantify the effects 
of blending. 

7. Simulations 

To better understand the effects of blending we 
have run extensive simulations of each of our NIC- 
MOS fields following the procedures of Stephens et 
al. (2001a). We create an artificial field to match 
each observed field, and measure it in exactly the 
same manner as the real frame. Since in the simu- 
lations we know both the measured and true mag- 
nitude of every star in the field, we can try to es- 
timate the true properties of the observed stellar 
population being modeled, free from observational 
effects. 

One of the goals of this work was to look for 
variations in the stellar populations with vary- 
ing galactocentric distance and bulge-to-disk ratio. 
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Table 5 
iC-BAND Luminosity Functions 



K 


Fl 


F2 


F3 


F4 


F5 


F170 


F174 


F177 


F280 


14.875 


0.00 


0.00 


0.00 


0.01 


0.00 


0.00 


0.00 


0.00 


0.00 


15.125 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


15.375 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


15.625 


0.01 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


15.875 


0.01 


0.00 


0.01 


0.00 


0.01 


0.00 


0.00 


0.00 


0.00 


16.125 


0.05 


0.00 


0.03 


0.01 


0.00 


0.00 


0.05 


0.02 


0.00 


16.375 


0.15 


0.01 


0.06 


0.00 


0.05 


0.03 


0.05 


0.09 


0.01 


16.625 


0.32 


0.01 


0.07 


0.05 


0.06 


0.06 


0.07 


0.26 


0.04 


16.875 


0.61 


0.02 


0.17 


0.09 


0.15 


0.11 


0.18 


0.50 


0.01 


17.125 


1.06 


0.02 


0.32 


0.08 


0.17 


0.14 


0.34 


0.53 


0.01 


17.375 


1.71 


0.00 


0.52 


0.21 


0.36 


0.19 


0.64 


1.25 


0.01 


17.625 


3.18 


0.06 


0.96 


0.50 


0.59 


0.60 


1.46 


1.79 


0.04 


17.875 


3.72 


0.07 


1.41 


0.69 


0.78 


0.73 


1.69 


2.48 


0.12 


18.125 


4.24 


0.03 


1.81 


0.83 


0.88 


0.64 


1.94 


2.81 


0.04 


18.375 


4.15 


0.15 


2.16 


0.98 


1.19 


0.88 


2.30 


3.13 


0.09 


18.625 


4.36 


0.18 


2.40 


1.15 


1.19 


1.38 


2.22 


3.18 


0.13 


18.875 


4.61 


0.19 


2.22 


1.16 


1.26 


1.12 


2.71 


3.18 


0.06 


19.125 


4.42 


0.10 


2.50 


1.41 


1.67 


1.18 


2.81 


3.66 


0.07 


19.375 


4.01 


0.14 


3.10 


1.50 


1.69 


1.43 


2.72 


3.49 


0.23 


19.625 


3.11 


0.25 


3.61 


2.04 


2.13 


1.71 


3.32 


3.22 


0.16 


19.875 


2.33 


0.28 


3.77 


1.95 


2.52 


2.08 


3.02 


2.91 


0.32 


20.125 


1.71 


0.35 


3.67 


2.61 


2.93 


2.49 


2.74 


1.98 


0.36 


20.375 


1.10 


0.38 


2.84 


2.55 


2.52 


2.28 


2.19 


1.01 


0.27 


20.625 


0.41 


0.37 


2.16 


2.47 


2.41 


2.53 


1.55 


0.60 


0.36 


20.875 


0.19 


0.39 


1.05 


1.79 


2.05 


1.85 


0.75 


0.26 


0.58 


21.125 


0.03 


0.33 


0.37 


1.01 


0.68 


1.00 


0.21 


0.07 


0.38 


21.375 


0.01 


0.17 


0.03 


0.17 


0.14 


0.30 


0.05 


0.00 


0.29 



However, the severe crowding, strong dependence 
of blending on surface brightness, and degeneracy 
between surface brightness and bulge-to-disk ratio 
make this question very difficult to answer. Thus 
one of the main purposes of our simulations is to 
determine whether all of our observations are con- 
sistent with a single stellar population. To make 
this determination, we have generated artificial 
frames using the stellar properties measured in 
some of the least crowded fields. If the measured 
differences between fields are just due to observa- 
tional effects, the simulations should exhibit the 
same differences. 

The simulations are complicated by the fact 
that the data come from four different instrument 
configurations, each with different exposure times, 
plate scales, and dither sizes. For each configu- 
ration, each dither starts as a blank frame hav- 
ing the appropriate noise characteristics. We then 
randomly add stars using the DAOPHOT ADDSTAR 
routine, until we have approximately matched the 
observed stellar density in the field being modeled. 
The PSFs used to add stars are the average of the 
PSFs determined from each field for each configu- 
ration, with any negative values in the model PSF 
set to zero. The ADDSTAR routine also incorpo- 



rates random Poisson noise into each star as it is 
added. 

The input stellar population was chosen to 
match the colors and luminosity function observed 
in the least-crowded bulge fields. The colors are 
the mean colors observed in fields 4 & 5, calcu- 
lated at 0.5 magnitude intervals. The input LF is 
a broken power law with a faint end slope of 0.278 
extending from — 5.7 < Mk < 5.0, and a bright- 
end slope of 1.100 from -7.4 < M K < -5.7. 
The faint-end slope was taken from the Galac- 
tic Bulge (DePoy et al. 1993), while the break- 
point and bright-end slope were determined from 
the NIC1 J-band LFs, since it was felt that they 
were the most robust against the effects of blend- 
ing. The artificial frames were then processed and 
measured in exactly the same manner as the real 
data, namely finding stars on a combined image 
with daofind, then measuring all dithers simul- 
taneously with ALLFRAME. 

The number of input stars was varied to ap- 
proximately match the number of detected stars 
on each frame, although the measured LF mor- 
phology was also taken into account for some of 
higher density fields. Table 6 lists the number of 
artificial stars input into each simulated frame, as 
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Fig. 10. — Comparison of M31 luminosity func- 
tions with that of Baade's Window (BW). The top 
panel shows the M31 disk (F2 & F280) LFs (solid) 
compared with the BW LF (beaded). The bottom 
panel shows the M31 bulge LFs (solid) overplotted 
on the BW LF (beaded). All LFs have been nor- 
malized in the range -7 < M K < -5.5. The BW 
LF is a composite of measurements made by Fro- 
gel & Whitford (1987) and DePoy et al. (1993). 
Note that the M31 disk LFs extend > 0.5 mag- 
nitudes brighter than the cutoff at Mk ~ —7.4 
seen in BW, while the M31 bulge LFs are in good 
agreement with the observations of BW. 

well as the number of stars recovered from both 
the real and simulated fields. 

The results of the simulations of a single cam- 
era (NIC2) and a single field (Fl) are summarized 
in Figure 12. This is the most crowded field, and 
hence the effects of blending are the most severe. 
The simulated J-band image is shown in the upper 
left, and is nearly indistinguishable from the real 
observations (Fig. 3). The resulting Mk,(J — K) 
CMD is shown in the upper right. The input 
stars form a narrow locus which blends into a 
line stretching (on this plot) from Mk — — 3 and 
(J - K) = 0.7 to M K = -7.4 and (J - K) = 1.48. 
The locus of measured stars is broader and shifted 
brighter and bluer compared to the input stars. 
As the simulations show, when the field is this 
crowded, none of our measurements are very ac- 
curate. 

The input and recovered luminosity functions 
are shown at the bottom of Figure 12. The J- 
band LFs on the left show that the measured LF 
is shifted by ~ 0.5 magnitudes brighter than the 
input LF. The LFs on the right show that, the K- 
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Fig. 11 — The number of M31 RGB stars within 
one magnitude of the RGB tip per resolution el- 
ement, based on the formulae of Renzini (1998). 
The left panel shows N(RGBT) per resolution el- 
ement as a function of the if -band surface bright- 
ness (magnitudes arcsec -2 ) for four different imag- 
ing resolutions. The right panel shows N(RGBT) 
per resolution element as a function of the distance 
from the center of M31 in arcminutes, based on 
the major-axis surface brightness measurements of 
Kent (1989). The transformation from Kent's r- 
band measurements to if -band surface brightness 
assume a constant color of (r — if ) = 2.9. 



band measurements are not as severely distorted 
as the undersampled NIC2 J-band observations. 

7.1. Simulation Results 

The first statistic we calculate is the differ- 
ence between the brightest star measured and the 
brightest star input into each artificial field. This 
gives a rough idea of the maximum amount of 
brightening one can expect in each field. Fig- 
ure 13 illustrates this difference for the J,H, and 
if -bands as a function of field surface brightness. 
This plot shows that there is little brightening due 
to blending in the artificial fields with lower sur- 
face brightnesses (fiK > 18). However, as the sur- 
face brightness increases, the amount of artificial 
brightening increases as well. 

The size of each circle in Figure 13 indicates 
whether the simulation is a NIC1 (small circles) 
or NIC2 (large circles) field. The higher resolu- 
tion and more finely sampled NIC1 observations 
are clearly less affected by blending, with only a 
few fields having deviations greater than 0.1 mag- 
nitudes. 

This figure also shows that the brightest K- 
band data are less affected by blending than the 
corresponding J-band measurements. For the 
most crowded fields (e.g. Fl, F177, F174, F3) the 
J-band brightening due to blending can be as high 
as 0.75 magnitudes, while in the if -band, blending 
is ~ 0.2 magnitudes less. This can also be seen in 
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the larger difference between the input and recov- 
ered LFs at J compared to K 7 and the blueward 
shift of stars in the CMD as seen in Figure 12. As 
mentioned above, we attribute the difference be- 
tween J & if to the undcrsampling in J and the 
bluer color of the underlying population. 

The delta-magnitude quantity previously calcu- 
lated is admittedly subject to small number statis- 
tics, since it is based on only a single star in each 
field. Ideally we would run each simulation a num- 
ber of times, and the average difference would be 
a much more robust estimator of the amount of 
brightening to expect in the real field. However, 
looking at the results of all the fields together, the 
direction and magnitude of the effect is clear. 




15 16 17 18 19 

Fig. 13. — The difference between the brightest 
measured and brightest input star in each simu- 
lated field plotted as a function of the field K- 
band surface brightness (mag/arcsec 2 ). The top 
panel shows the J-band difference for the NIC1 
and NIC2 fields, the bottom two panels show the 
H- and if-band difference for the NIC2 fields. The 
Field ID is indicated in the middle of each cir- 
cle, where the large circles denote the larger NIC 2 
fields, and the smaller circles the smaller NIC1 
fields (J-band). 



Another interesting, and hopefully more robust, 
quantity to calculate for each field is the ratio of 
the number of "bright" stars measured compared 
to the number of stars input to the same bright- 
ness. This quantity is very important, for example 
when using AGB stars to assess recent star for- 
mation. We plot this ratio as a function of the 
field surface brightness in Figure 14. We chose 
Mj < -5, M H < -6, and M K < -6.5 as the cri- 
terion for a star to be considered "bright" . These 
limits are fairly arbitrary, however if chosen to be 
much brighter, then some of the fields will have no 
stars input that bright, and if much fainter, some 
fields will not be complete to that level, and we 
will be measuring completeness instead of blend- 
ing. 

Figure 14 shows the ratio of measured to input 
bright stars for all of our simulated frames. The 
large circles represent NIC2 fields and the small 
circles NIC1 fields. The lowest surface brightness 
fields (F2 & F280) have nearly equal numbers of 
measured and input stars, i.e. minimal blend- 
ing. However, as the surface brightness increases, 
the points begin to move up off of the dashed 
line indicating a ratio of unity. This upturn is 
a function of wavelength, but in general occurs at 
fj,K ~ 17 mag/arcsec 2 for the NIC2 observations, 
and [Ik ~ 15.5 mag/arcsec 2 for the NIC1 obser- 
vations. In the simulations of the brightest fields 
we measure about twice as many bright stars in 
NIC2, and 1.25 times as many in NIC1 compared 
to the number which were actually input into the 
simulation. 

Of course the ratio of the number of measured 
to input stars is dependent on exactly where one 
draws cutoff magnitude. At brighter cutoff mag- 
nitudes the ratio goes to infinity when there are 
no stars input as bright as we measure. Choosing 
a fainter cutoff both dilutes the number of blends, 
and causes faint blends to be lost due to incom- 
pleteness. 

In summary, both Figures 13 and 14 show sim- 
ilar structure, with a sharp increase in blending 
between [Ik ~ 16 — 17. Obviously the amount of 
blending one can withstand depends on the scien- 
tific goals, however one must be very careful inter- 
preting results from such data. We are even skep- 
tical of some of our own measurements of stars just 
above the tip of the AGB, which only occur in the 
most crowded fields where the surface brightness 
is greater than fiK ~ 16 magnitudes arcsecond~ 2 
(however see §9). 
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Fig. 12. — Simulation of the NIC2 field Fl. The combined J-band image is illustrated in the upper left. 
The CMD showing both input (narrow locus of stars) and measured stars is in the upper right. The J and 
if-band input (solid) and measured (beaded) luminosity functions are shown in the bottom two panels. 
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Fig. 14. — The ratio of the number of bright stars 
measured in each simulated field to the number 
of bright stars input plotted as a function of field 
surface brightness (mag/arcsec 2 ). The top panel 
shows this ratio for all the J-band frames, where 
we have chosen Mj < — 5 as the criterion to be 
counted as a bright star. In this plot the NIC1 
fields (small circles) are less affected by blending 
than their NIC2 counterparts (large circles). The 
middle panel shows the ratio of measured to input 
bright stars in all of our iJ-band frames, where we 
only count stars with Mh < —6, and the bottom 
panel gives the ratio for the if -band frames, using 
M K < -6.5. 



8. Comparison with Previous Observa- 
tions 

8.1. Rich, Mould & Graham (1993) 

The groundbased observations of Rich, Mould 
& Graham (1993, hereafter RMG93) are the foun- 
dation of the current work. Theirs was the first 
study in M31 which systematically attempted to 
measure the stellar properties over a range of 
galactocentric distances and bulge-to-disk ratios. 
For that reason our NIC2 fields were chosen to be 
the same as those in RMG93. RMG93 took their 



observations on 1992 Aug 30 - Sept 1 with the 
Palomar IR imager on the Hale 5m telescope. This 
instrument was outfitted with a 58 x 62 pixel InSb 
detector with 0.313" pixels. Each of their 18" x 19" 
fields were observed for 75 seconds through both J 
and K filters. Using offsets equal to half the field 
of view, they obtained 25 frames that were later 
assembled into 72.6" x 77.6" mosaics. Thus the 
central 60" of their mosaics have total integration 
times of 150s in each filter. The seeing is ~ 1" on 
these mosaics. 

The central - 60" of Field 1 from RMG93 
is shown on the left side of Figure 15. In or- 
der to match up our corresponding NIC2 image 
shown in the right side of Figure 15, we first re- 
binned our image by a factor of 4.1 to go from our 
0.0757" pixels to the groundbased 0.303" pixels, 
then smoothed the rebinned image with a Gaus- 
sian kernel to match the 1" seeing in the ground- 
based images. This rebinned and smoothed inter- 
mediary NICMOS image is shown in the center of 
Fig 15. 

In order to better understand the relationship 
between the groundbased photometry and our 
NICMOS measurements we have performed a star- 
by-star comparison of the objects measured by 
RMG93 in Field 1. As is obvious from Figure 15, 
none of the objects seen from the ground corre- 
spond to single stars in the NICMOS image. How- 
ever, if we simply take our brightest measured star 
nearest the RMG93 object as the center of the 
clump which composes their object, we can study 
the composition of that clump. 

As an example consider RMG93 star 95. This 
object is located just above and right of the center 
of the groundbased image. RMG93 measure this 
object as K = 15.15, however the star we match 
this object with has (the first entry in Table 7) 
K = 16.63. The next star has K = 16.35 and 
lies 0.15" away. If we include all stars within a ra- 
dius of 0.5" we should get approximately the same 
amount of flux as measured by RMG93 viewed 
through 1" seeing. Table 7 lists the NICMOS mea- 
sured stars and their radius from the assumed cen- 
ter of the groundbased clump (columns 2 & 3) and 
the magnitude of the running sum of their flux 
(column 4). Going down through this table, we 
eventually add enough stars to reach the ground- 
based measurement of K = 15.15 at r ~ 0.52". Of 
course this radius of equal measurements varies 
from star to star, depending on how PSFs were 
fit to the groundbased blends. Using 52 "stars" 
matched with the RMG93 observations, we find 
this average radius to be 0.47" in K and 0.35" 
in J (see Appendix A). Of the 13 stars with 
matches in field 1, the average difference between 
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Table 6 
Numbers of stars 



Recovered 



Field 


Input 


Real 


S i mu 1 a.t ion 


NIC1 








Fl 


900000 


4037 


3277 


F2 


25000 


348 


521 


F3 


400000 


3102 


3088 


F4 


150000 


2208 


2411 


F5 


250000 


2323 


2405 


F170 


300000 


2154 


2055 


F174 


950000 


4847 


4970 


F177 


950000 


4861 


4970 


F280 


50000 


448 


619 


NIC2 








Fl 


3000000 


4037 


3594 


F2 


80000 


314 


597 


F3 


1500000 


3128 


3498 


F4 


700000 


2066 


2598 


F5 


800000 


2257 


2555 


F170 


900000 


1986 


2496 


F174 


2000000 


2923 


3156 


F177 


4000000 


3221 


3010 


F280 


110000 


256 


487 




Fig. 15. — Comparison between groundbased if-band image of Fieldl (left) and the NIC2 F222M image 
(right). The center image is a rebinned and smoothed version of the NIC2 image used to help match up the 
observations. 
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the RMG93 measurement and the brightest NIC- 
MOS star measured within 0.47" is —1.41 magni- 
tudes = 0.30). 

Figure 16 shows a comparison between the NIC- 
MOS measured LFs (solid lines) with RMG93's 
groundbased LFs (beaded lines). The RMG93 
LFs are the measured numbers of stars matched in 
their J and X-band images placed into 0.25 mag- 
nitude bins (column 2 of their Table 8). We have 
normalized these LFs by multiplying by 0.063 to 
compensate for the larger area of the groundbased 
images (5634 arcsec 2 ) compared to the NIC2 area 
(355 arcsec 2 ). If the normalized RMG93 LF falls 
below a value of one (dashed line) it can be taken 
as the probability that NICMOS would find a star 
that bright if they exist. 

Looking at the comparison in Figure 16 we see 
that the amount of disagreement is strongly corre- 
lated to the field surface brightness, which is listed 
in the upper left of each panel. Field 1 exhibits 
the worst case of blending, where RMG93 finds a 
significant number of stars a magnitude brighter 
than we measure. Field 3 is not as severe, how- 
ever the groundbased observations predict that we 
should find many more bright stars up to ~ 0.5 
mag brighter than we see. The three lowest sur- 
face brightness fields (F2, F4, F5) are roughly con- 
sistent with the NICMOS measurements. In every 
field the RMG93 LF trails off to very bright mag- 
nitudes, and even though there is a very low prob- 
ability of finding such bright stars in one of our 
small NIC 2 fields, it seems clear that these are 
most likely just severe cases of blending. 

8.2. Rich & Mould (1991) 

Rich & Mould (1991, hereafter RM91) obtained 
the first infrared color- magnitude diagram of an 
M31 bulge field, 3.65' from the nucleus and along 
the major axis. They observed J and K-band mo- 
saics of nine non-contiguous 18" x 19" fields using 
the Palomar IR imager on the Hale 5m telescope, 
yielding a total area of ~ 3110 square arcseconds. 

The RM91 field lies close to our F3 field, and 
we compare the luminosity functions by taking the 
RM91 counts from column 4 (CMD) of their Table 
3. We then normalize the RM91 LF by multiplying 
by the ratio of the NICMOS field F3 area (355 
arcsec 2 ) to the RM91 field area. The resulting 
comparison is shown in Figure 17. 

Figure 17 shows remarkable agreement between 
the bright ends of the NICMOS and RM91 lumi- 
nosity functions. Had the M31 bulge population 
been normalized to Baade's Window in RM91, the 
apparent extended giant branch would have been 
far less prominent. However an accurate normal- 




Fig. 16. — Comparison between our NIC2 mea- 
sured LFs (solid lines) and the LFs measured 
by Rich, Mould & Graham (1993) (beaded lines) 
scaled to match the NICMOS field area. The K- 
band surface brightness of each field is listed under 
the field label in the upper left of each panel. 

ization is difficult without the faint end complete- 
ness provided by the NICMOS images. 

Figure 17 also shows a discrepancy between 
the luminosity functions of RM91 and RMG93, 
where the earlier RM91 groundbased data seems 
to match the bright end of the NICMOS LF, while 
the later RMG93 data appears much more affected 
by blending. 

One possible explanation for this difference 
seems to lie in the analysis of the data. While 
both datasets were acquired with the same tele- 
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Table 7 
RMG93 Star 95 (K = 15.15) 



N 


Distance* 


K 


Sum 


1 


0.01 


16.63 


16.63 


2 


0.15 


16.35 


15.73 


3 


0.25 


17.87 


15.59 


4 


0.25 


17.24 


15.37 


5 


0.44 


18.68 


15.32 


6 


0.47 


18.49 


15.26 


7 


0.48 


18.75 


15.22 


8 


0.51 


19.47 


15.20 


9 


0.52 


18.84 


15.16 


10 


0.52 


19.27 


15.14 


11 


0.57 


19.47 


15.12 


12 


0.59 


17.09 


14.96 



a Distance in arcseconds. 



scope and instrument, the RM91 data were ana- 
lyzed on a frame-by-frame basis, while the RMG93 
data were assembled into a large mosaic before 
analysis. Therefore the RM91 data retained its 
original image quality, while the RMG93 image 
quality, due to difficulties in perfectly registering 
the frames, was reduced to perhaps even worse 
than that of the worst image. 

In light of our simulations, it is still difficult to 
understand the apparent agreement between the 
bright ends of the RM91 and NICMOS luminosity 
functions. Our simulations (for their quoted 0.6" 
seeing) would still lead to the inescapable conclu- 
sion that the measured magnitudes of the bright- 
est stars in RM91 must still suffer from crowding. 
It is also possible that the images may have been 
so undcrsampled (0.31" pixels) that that the ac- 
tual seeing was better than quoted, though RM91 
state that 0.6 arcsec seeing was reached in only 
one of the images. As the original frames are not 
available, we are not able double check the mea- 
surements or to explore the issue in further detail. 

8.3. Davidge (2001) 

Davidge (2001) has recently obtained JHKs 
images of a bulge field 2.6' SW of the nucleus of 
M31 using the 3.6 meter Canada-France-Hawaii 
Telescope (CFHT). With the help of adaptive op- 
tics, his images achieve a FWHM of 0.35". His 
photometric uncertainties include 0.05 magnitudes 
in the aperture correction and 0.03 magnitudes in 
the zero-point. The surface brightness of his field 
(00 h 42 m 45.1 s , +41°13'31.3", J2000) is fj, K ~ 15.5 
mag/arcsec 2 based on the measurements of Kent 
(1989) and assuming (r - K) = 2.9. 




K 

Fig. 17. — Comparison between our F3 if -band 
luminosity function (solid line) and the LF mea- 
sured by Rich & Mould (1991) (beaded line). Also 
plotted is the field 3 measurement of Rich, Mould 
& Graham (1993). We have scaled the ground- 
based data to match the NICMOS field area. 

Figure 18 shows a comparison between the LFs 
of the Davidge (2001) field and the two NICMOS 
fields with bracketing surface brightnesses: F174 
and F177, which have \ik — 15.8 and [Ik — 15.4 
mag/arcsec 2 respectively. The Davidge LF should 
lie between the two NICMOS LFs, however the 
Davidge LF is instead shifted ~ 0.5 magnitudes 
brighter. 

Davidge has suggested that the difference be- 
tween his measurements and our HST-NICMOS 
observations (Stephens et al. 2001b) is due pri- 
marily to calibration. In support of this conclu- 
sion he cites simulations which indicate that the 
effects of blending on his observations are at most 
0.1 mag over what would be measured with the 
NIC2 resolution. However, these simulations used 
simple Gaussian PSFs and included only ~ 10000 
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stars (< 10 stars/arcsec 2 ). They show only a few 
severely blended stars which he claims are "eas- 
ily identifiable" . In reality blending is a stochastic 
phenomenon, involving millions of stars, and pro- 
ducing a continuum of blending which can be very 
difficult to detect and quantify 

As shown by Figure 11, making observations at 
/j,k = 15-5 mag/arcsec 2 at Davidge's 0.35" reso- 
lution will give, on average, ~ 0.4 NRGBT stars 
per resolution element. Thus in any given reso- 
lution clement there is a (0.4) 2 = 0.16 probabil- 
ity of having 2 NRGBT stars, compared to the 
(0.07) 2 = 0.049 probability with NICMOS. Thus 
the most likely cause for the difference between 
the Davidge (2001) LF and our NICMOS LFs in 
Figure 18 is indeed due to blending. 

As wc discussed in Section 7, we are certainly 
not claiming the comparison of our observations 
with Davidge (2001) is a case of right and wrong; 
but rather a case of wrong and wrong. Both sets 
of observations are affected by blending, and the 
stars just above the AGB tip in our most crowded 
fields may in fact be artificially brightened by sev- 
eral tenths of a magnitude. 
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Fig. 18. — Comparison between the M31 bulge 
LFs measured by NICMOS: F174 (solid line) 
and F177 (dashed line), with that measured by 
Davidge (2001) (beaded line). The surface bright- 
nesses of each field is: 15.8, 15.4, and 15.5 K mag- 
nitudes per arcsec 2 for F174, F177, and Davidge's 
field respectively. 

9. Bright Stars 

We have cautioned that some of the bright stars 
in our fields may be blends of fainter stars, however 



there exists a population of bright stars which are 
real. Inspection of the images shows that these 
stars are obvious point sources, and occur over 
the entire range of M31 stellar densities. These 
stars are some of the brightest and bluest which 
we have observed, with 16.68 > K > 13.75 and 
0.45 < ( J - K ) < 0.75, and arc most likely fore- 
ground Milky Way stars. Here we provide a brief 
discussion of their properties to ensure that they 
are not confused with a population of young M31 
bulge stars. 

The CMD of brightest stars measured in all of 
the NIC2 frames is shown in Figure 19. Here the 
AGB is located between 1 < (J — K) < 2 and 
extends up to the curved line at M^oi = — 5. Stars 
directly above the AGB are most likely blends of 
fainter stars. However the bright stars bluer than 
(J — K) = 0.8 are indeed real stars. If they were at 
the distance of M31 they would all have bolomctric 
magnitudes brighter than M^i = —5.5. 

We find a total of 8 foreground stars in the NIC2 
fields: 2 in F2, 1 in F3, 2 in F4, 1 in F5, 1 in F170, 
and 1 in the Gl field, which has not been otherwise 
analyzed in this paper because of the lack of non- 
cluster stars. We find none in the NIC2 Fl, F174, 
F177, or F280 fields. 

Searching the NIC1 observations, we use the 
fact that the faintest foreground star in the NIC2 
group has J = 17.43. Thus if we assume that any 
star with a NIC1 J-band magnitude brighter than 
J = 17.5 is also a MW contaminant, wc find 2 in 
the NIC1 F170 field, and 1 in the F177 field. None 
of the other 8 NIC1 fields (including Gl) have any 
stars this bright. These three stars are also clearly 
separated from the rest of the NIC1 AGB LF by 
a > 0.25 magnitude gap. 

The radial distribution of the bright foreground 
stars is illustrated in Figure 20. Here wc plot the 
surface density measured in each field as a function 
of radial distance from the center of M31. The up- 
per and lower limits are one-sigma confidence in- 
tervals, calculated using the small-number approx- 
imation formulae of Gehrels (1986). The distribu- 
tion shows no trend with radius, and is certainly 
not correlated with either of the steeply dropping 
surface brightness profiles of M31's bulge or disk, 
illustrated by the dotted and short-dashed lines 
respectively. 

In total, we find 8 bright foreground stars in 
the 10 NIC2 fields. These fields have a combined 
area of 0.95 arcmin 2 , giving a surface density of 
8.4 arcmin 2 . Although we don't have color infor- 
mation in the NIC1 frames, there are three very 
bright stars with J < 17.5. The total area of the 
NIC1 fields is 0.56 arcmin 2 , giving a surface den- 
sity of 5.4 stars arcmin -2 for these bright stars. If 
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Fig. 19. — The top end of the combined CMD 
for all (NIC2) frames. Note the clear separation 
between the red M31 AGB stars and the bluer 
foreground MW stars. The curved line illustrates 
Mb i = —5, thus «/these stars were at the distance 
of M31, they would all have bolometric magni- 
tudes over —5.5. 

we assume that all 11 stars are of the same pop- 
ulation, the average surface density is 7.3 stars 
arcmin -2 , which is shown by a long-dashed line in 
Figure 20. 

A comparison with the estimated number of 
field stars by Ratnatunga & Bahcall (1985) re- 
veals that our measurement of 7.3 stars arcmin -2 
is extremely high. The measured mean color of 
(J — K) = 0.6 corresponds to (V - K) ~ 2.2 for 
either dwarfs or giants, which means that these 
stars most likely have 18.88 > V > 15.95. How- 
ever, Ratnatunga & Bahcall (1985) predict ~ 0.8 



field stars / arcmin 2 between V = 15 and V = 19 
toward M31 . Given the combined area of our NIC1 
and NIC2 images (1.5 arcmin 2 ), we should have 
found approximately one field star, not 11. 

Unfortunately there have been few surveys of 
bright field stars toward M31 to verify the model 
predictions. Ferguson et al. (2002) recently per- 
formed a large scale survey of M31 looking for 
substructure in the halo and disk of M31. In- 
tegrating over all their magnitudes, 18 <i< 23, 
which correspond to roughly 19 < V < 24, they 
estimate 13,000 - 20,000 Galactic foreground stars 
dcg~ 2 or between 3.6 - 5.56 stars arcmin -2 . This 
is very close to the prediction of Ratnatunga & 
Bahcall (1985) for their magnitude range (~ 3.5 
/arcmin 2 ), but the stars we have observed mostly 
have V < 19, and thus are too bright to be in- 
cluded in their survey. 

Thus while the crowded fields have many blends 
at Mb i ~ —5, it is clear that the brightest and 
bluest stars, with K < 16.8 and (J - K) < 0.8, 
are real. We find 11 of these stars which is over 
a factor of 10 greater than what is predicted by 
Galactic models. However it seems very unlikely 
that these stars are associated with M31 as their 
surface density does not scale with the surface 
brightness of M31, and we find one as far out as 
the globular cluster Gl, 34 kpc from the center of 
M31. 

10. Discussion & Conclusions 

We have analyzed the stellar populations of 
M31 using 9 sets of adjacent HST NIC1 and NIC2 
fields, with distances ranging from 2' to 20' from 
the nucleus. These observations are the highest 
resolution IR measurements to date, and provide 
some of the tightest constraints on the maximum 
luminosities of stars in the bulge of M31. 

Analytic estimates of the effects of blending on 
our observations, where we calculate the number 
of RGB stars within 1 magnitude of the RGB tip 
as a function of position in M31, indicate that sim- 
ulations are required to accurately interpret our 
observations. We thus perform extensive simula- 
tions of each of our NICMOS fields following the 
procedures of Stephens et al. (2001a). These sim- 
ulations show that for the most crowded fields we 
can expect the brightening due to blending to be 
as high as 0.75 magnitudes in J and 0.55 mag- 
nitudes in K. They also show that the ratio of 
measured to input bright stars is a strong func- 
tion of the field surface brightness. In the high- 
est surface brightness field we measure about 25% 
more bright stars than were input with NIC1, and 
about twice as many with NIC2. 
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Fig. 20.— The number of bright foreground stars 
per square arcminute as a function of radial dis- 
tance from the nucleus of M31 (left axis). The 
long dashed line at 7.3 shows the mean taken over 
all fields (11 stars over 1.5 arcmin 2 ). We have 
also overplotted the major axis surface brightness 
profiles of the bulge (dotted line) and disk (short 
dashed line) of M31 as measured by Kent (1989) 
(right axis). Here the surface brightness offset is 
arbitrary, but the scale is set to match the number 
counts, so that if the bright stars were associated 
with a population in M31, they should follow the 
radial surface brightness profile of that population. 



in his lower resolution observations, rather than a 
calibration error as claimed by Davidge. 

We also find an unusually high number of bright 
blueish stars in our fields. In all 20 (NIC1 + NIC2) 
fields we find 11 stars which are uncorrelated with 
the surface brightness distribution of M31, and ap- 
pear to be foreground Milky Way stars. However, 
the implied surface density of 7.3 arcmin^ 2 is over 
a factor of 10 higher than is predicted by Galactic 
models. 

Support for this work was provided by NASA 
through grants GO-7826 and GO-7876 from the 
Space Telescope Science Institute. Thanks to Pe- 
ter Stetson for supplying, and helping us with his 
ALLFRAME photometry package. JAF thanks Dr. 
Sean Solomon for providing Visiting Investigator 
privileges at DTM/CIW. Also, many thanks to 
Sergio Ortolani for giving a very helpful referee's 
report. 



All of the bulge luminosity functions are consis- 
tent with a single, uniform bulge population. This 
is based on the observation that the small differ- 
ences we see are correlated with surface bright- 
ness, and that the simulations predict similar dif- 
ferences. We note however, that our simulations 
only verify that a single LF combined with blend- 
ing can produce the observed field-to-field differ- 
ences. Thus our observations are consistent with 
a single bulge LF, but without higher resolution 
data, small field to field differences cannot be ruled 
out. 

The tip of the RGB in M31 is clearly visible at 
Mboi ~ —3.8, and the tip of the bulge AGB ex- 
tends to Mk ~ —8. This AGB peak luminosity 
is significantly fainter than previously claimed. A 
comparison with the measurements of Rich, Mould 
& Graham (1993), which guided our choice of the 5 
pointed observations, indicates that their brightest 
stars are most likely severe cases of blending. In 
a comparison of our F174 and F177 LFs with the 
recently observed bulge field of Davidge (2001), 
we find that the ~ 0.5 magnitude difference be- 
tween his and our LFs is due entirely to blending 
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A. NIC1 Transformation 



In order to compare observations made with the two different NICMOS cameras with each other and with 
groundbased observations wc must first convert all measurements to a common photometric system. The 
transformation of NIC2 to the groundbased CIT/CTIO system has already been calculated and published 
(Stephens et al. 2000). However, NIC1 lacks a formal transformation to any groundbased system. 

As a first attempt to transform NIC1 to a groundbased photometric system, we applied the NIC2 trans- 
formation from Stephens et al. (2000). Using their calibration keywords, and assuming (J — K) = 1 for the 
color term, we applied the corresponding offset to all of our NIC1 photometry. However, this transformation 
yielded large discrepancies (up to ~ 0.5 magnitudes) between the luminosity functions measured in corre- 
sponding NIC1 and NIC2 field pairs, such that the NIC1 fields appeared too faint. A comparison between 
the STScI calibrated NIC1 and NIC2 F110W luminosity functions shows that the lower surface brightness 
fields should show nearly perfect agreement, while the more crowded fields should show discrepancies of no 
more than ^0.25 magnitudes. 

As an alternative method to transform NIC1 to a groundbased system, we considered the observations of 
Rich, Mould & Graham (1993). Our NIC2 observations were chosen to be centered on the RMG93 fields, 
and since the NICI and NIC2 focal planes are so close together (17.5" between field edges) and NICMOS 
was rotated ~ 45 degrees from North when we took most of our images, NICI falls on the lower left (SE) 
corner of the RMG93 images (see Figure 2). By rebinning and smoothing our NICI images we were able 
to match up 8 "stars" with RMG93. Of course these aren't really single stars, but rather clumps of many 
stars. However, by estimating how many stars are in the clumps measured by RMG93, we were able to use 
their observations to transform ours to the CIT/CTIO system. 

Going back to our NIC2 observations, whose calibration we trust, we matched up 52 "stars" with RMG93. 
Using the J-band observations, we determined that we can make NICMOS agree (±0.3 magnitudes) with 
RMG93's measurements if we sum up all the NICMOS measured stars within a 0.35" radius around what 
we estimate to be the center of the RMG93's "stars". 

To determine the NICI transformation, we first calibrate our photometry using the most recent header 
keywords listed in the NICMOS Data Handbook v5.0 (Dickinson et al. 2002) namely: photfnu = 2.358E-6 
Jy sec / DN, and fnuvega = 1773.7 Jy. We then sum the flux of all NICMOS stars measured within 0.35" 
of the RMG93 ccntroids. The resulting difference between the NICI magnitudes and RMG93 is 0.36 ± 0.17 
magnitudes using all 8 "stars", or 0.42 ± 0.09 magnitudes using a sigma-rejected sample of 7. 

In summary, we apply a —0.42 magnitude offset to our NICI F110W magnitudes to approximately trans- 
form them to the J-band of the groundbased CIT/CTIO system. This is in contrast to the NIC2 photometry 
which was transformed using the equations of Stephens et al. (2000), which includes a color term in each 
band. 
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